Influence of Disorder on Electron-Hole Pair Condensation in Graphene Bilayers 
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Graphene bilayers can condense into a state with spontaneous interlayer phase coherence that 
supports dissipationless counterflow supercurrents. Here we address the influence of disorder on the 
graphene bilayer mean-field and Kosterlitz-Thouless critical temperatures and report on a simple 
criteria for the survival of pair condensation. 
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Introduction — Superconductivity is a spectacular and 
potentially useful condensed matter phenomena. Unfor- 
tunately its occurrence has so far been limited to rel- 
atively low temperatures, even in the case of high-T c 
materials. The effective attractive interaction between 
two electrons that is imperative for superconductivity re- 
lies on some artifice that diminishes the role of the bare 
strongly repulsive Coulomb interaction. It has recently 
been suggested that higher temperature superfluid be- 
havior can be realized in graphene bilayers when con- 
ditions are favorable for Coulomb-driven electron-hole 
pair formation J3, Q , and that interesting new superfluid 
phenomena^, 0] will occur in the resulting counterflow 
superfluids. One obstacle which stands in the way of re- 
alizing this state is the inevitable presence of unintended 
disorder. In this Letter we derive a simple criterion for 
the survival of spontaneous interlayer coherence and con- 
terflow superfluidity in graphene bilayer systems which 
is informed by current understanding [5| of the disorder 
present in these truly atomic Q two-dimensional electron 
systems. 

The system we consider [3, Q consists of two parallel 
graphene layers embedded in a dielectric. The thin bar- 
rier separating the layers suppresses inter-layer tunnel- 
ing jjjj]. External gates can shift the Dirac cones and 
transfer charge between the layers. Spontaneous inter- 
layer coherence is most likely to occur [3, 0] when the 
energetic Dirac-cone shifts are equal and opposite so 
there is perfect nesting between the electron Fermi sur- 
face of the high-density layer and the hole Fermi surface 
of the low-density layer, i.e. when the bilayer is neu- 
tral. Nearly perfect nesting in neutral bilayers is guar- 
anteed by the nearly-perfect particle-hole symmetry of 
graphene's 7r-bands. At sufficiently low temperatures the 
system is driven to condensation by the Cooper instabil- 
ity. Graphene is the ideal system for bilayer electron-hole 
condensation Q because its bands are particle-hole sym- 
metric, because the lack of a gap assists charge transfer, 
and because its linear dispersion and truly-2D character 
increase electron- hole pairing energy scales at a given car- 
rier density. The stiffness of the interlayer phase in the 
condensed state facilitates counterflow superfluidity and 
in separately contacted bilayers supports novel transport 



phenomena which are still relatively unexplored 0, 0]. 

In this work we study the influence of disorder on pair- 
condensation in graphene bilayers. We find that disorder 
suppresses both the mean-field-theory critical tempera- 
ture T c and the KT temperature T KT , but that sponta- 
neous coherence survives when 



< 



< 
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Here is the density of the Coulomb scatterers that 
dominate 0, Q the disorder present in graphene sheets 
on dielectric substrates, n — kp/n is the external electric 
field induced carrier density in each layer, k F is the Fermi 
momentum, and d is the separation between graphene 
layers. The right inequality in Eq.([T]) expresses the re- 
quirement that interlayer interactions be comparable to 
intralayer interactions discussed in earlier work[l| while 
the left inequality places a limit on the allowable disorder 
strength. A window in carrier density for spontaneous co- 
herence exists provided that mcPir < 1. Since is typ- 
ically in the range between 10 11 cm -2 and 10 12 cm~ 2 for 
graphene on Si02 substrates, layer separations less than 
about lOnm are required for spontaneous interlayer phase 
coherence. In the body of this Letter we explain the de- 
pendence of T c and T KT on the strength of disorder and 
derive condition JTJ). 

Mean-Field- Theory with Disorder — In order to focus on 
the role of disorder in spontaneous coherence we sim- 
plify the mean-field-theory by neglecting the inessential 
role of the full valence band of the high-density layer 
and the empty conduction band of the low-density layer. 
Furthermore we incorporate the effects of the intra-layer 
Coulomb interaction only through an implicitly renor- 
malized Fermi energy e F . In the same spirit, the momen- 
tum dependence of the inter-layer interaction is replaced 
by an energy cutoff v/d (K = 1) where v is the velocity 
of the Dirac quasi-particles. The influence of disorder is 
incorporated here using the self-consistent Born approx- 
imation (SCBA). 

The basic building block in this analysis is the retarded 
Green function matrix G r aa , which satisfies the Dyson 
equation 
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In Eq.([2]) £fc = vp(k — kp) is the band energy, and the 
Pauli matrices act on layer labels aa' . The disorder and 
pairing self energies are respectively 

Sw(k,u) = n^lS^Re {K_ prT U^ pa ,} G^,(p,u;); 



and 
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where U qa is the disorder potential in layer a and 5k P = 
(l + e 4 ^ k_0p ^) /2 is graphcne's chiral form factor. The 
overbar in Eq.(4) denotes the average over disorder and 
g stands for the inter-layer interaction coupling constant. 
Since elastic impurity scattering occurs in the vicinity of 
the Fermi surface we can assume that U = U(0), i. e. that 
the scattering rate depends only on the angle between the 
incoming and outgoing momenta. In a similar spirit we 
also approximate the density-of-states by its value at the 
Fermi energy to obtain 



At x 



A 2 



where 



(w,A) = (w,A) + 



2r s ' 2t d J y/T~^' 



u = u/A, and 
1 
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(5) 



(6) 



(7) 



We have assumed here that the two layers are identical 
and set tu = T22 — t s (S = same layer) and T12 = T11 = 
t d (D=different layer). 

We now estimate t s and r D . The measured transport 
properties of graphene layers strongly suggest 
that disorder is dominated by Coulomb scatterers[5, 
9] near the interfaces between the dielectric and the 
graphene sheets. We therefore expect inter-layer disor- 
der potential correlations to play an inessential role and 
set t d — > 00. (Including a finite value of r D would not 
complicate our analysis in any way.) After accounting for 
the difference between transport and scattering lifetimes: 
Ttr/rs ~ 2 [n| we follow Ref.[0| in estimating that 
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We use this equation below to express r s in terms of n$ . 
Coherence Criterion — There is a complete analogy be- 
tween the Green function ((5]) and that of an electron in a 
superconductor with magnetic impurities once the pair- 
breaking parameter is identified with S = l/2r s + l/2r D . 



The various ways in which disorder influences coher- 
ence in bilayers may therefore be understood by bor- 
rowing results from Abrikosov-Gorkov theory [ll], [lij ]. 
In the condensed state the density of states (DOS) of 
a disordered bilayer graphene is v(u>) = ^Imu where 
( = S/A(T) and v = kp/irv is the DOS of the normal 
system. Disorder smoothes the square root singularity of 
a clean system. Furthermore it reduces the energy gap 
to A(T)[1 - ( 2 /3]3/2 for £ < 1. For C > 1 the spectrum 
is gapless. The mean field critical temperature T c in the 
presence of disorder is given by 11 1 



log 



<::0 



= H{a) 



(9) 



where TC(a) = tp (i + —J — ip (4) with ip being the 
di-Gamma function, a = 5/2T c , and T c q is the mean 
field critical temperature of the clean system. For weak 
disorder T c — T c q — ttS/A, i.e. the critical temper- 
ature decreases linearly with m. For strong disorder 
T 2 = 6 {S/irf log (7rT c0 /2£r) where T = 1.76. When 
5 > 0.8T c o the spectrum becomes gapless, however inter- 
layer phase coherence continues to exist. Eventually dis- 
order destroys superfluidity altogether when 5 > 0.88T n n- 
This last result of Abrikosov-Gorkov theory 



11] 



can 



be transformed into a more useful form by noting that 
for graphene bilayers in the strongly interacting regime 

(*B = 1) 



Too 



0.1e F 
kpd 



(10) 



Combining Eq.([10]) and Eq.([8]) yields the left inequality 
in (Q}. 

It is difficult to estimate T c q in the weakly interact- 
ing regime when k F d ^> 1. In that regime screening and 
other induced interaction effects are expected to signifi- 
cantly reduce the critical temperature. Precise estimates 
lie beyond the scope of mean field theory and present 
an interesting theoretical challenge. Nevertheless using 
Eq. ([9]) we can find that coherence between the two layers 
exists as long as 



> 
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For weak disorder, a <C 1, condition (fTT|) is equivalent 

to n > (jj^f-) e _7rQ . In the opposite limit (a>l) the 

system will condense for n > (ffj) 2 e~ 7I ' 2 / 12Q2 . 
Kosterlitz- Thouless Temperature — In two dimensions su- 
perfluidity is destroyed at the KT temperature by vor- 
tex anti-vortex proliferation. To estimate T KT we use 
the Kosterlitz- Thouless equation p s (Tkt) — {^/^)Tkt 
with p s (T) in the presence of disorder calculated from 
Abrikosov-Gorkov theory. The phase stiffness is deter- 
mined from the counterflow current j D generated by a 
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uniform gradient of the relative phase between the two 
layers and from the relation 2p s A = j D (A) — j D (A = 0). 
The phase gradient 2A is obtained by perturbing the two 
layers with a pair of constant vector potentials that are 
equal in magnitude but have opposite signs. The sub- 
traction of j D (A = 0) in the last relation is required[lj 
by a pathology of the Dirac model which implies that j D 
does not vanish, as required by gauge invariance, in the 
normal (A = 0) state. We find that 



Ps{T) = Pi 0) -^E / ^n F (e)ImTr[G r (k,cu)] 2 

= P M[i+ Jdu,^ J de^,e,A,o],(i2) 

where p s (0) = e F /47r is the zero temperature phase stiff- 
ness of a clean system fli and 



^(w,£,A,C) = --Im / 
* J-< 



e 
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It is essential that the integration over lj' in Eq. I|12p pre- 
cedes the integration over £. Vertex corrections to p s (T) 
vanish when the disorder potential U = U(6) is taken to 
be independent of the magnitudes of the incoming and 
outgoing momenta. 

In the clean limit J 7 is a sum of two delta functions 
centered at u = ±\Ji 2 + A 2 . The phase stiffness is then 



(0) 



du>— —v(uj). 
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The second term in (|14p captures the reduction in p s 
at finite temperature due to the thermal excitation of 
quasi-particles. It can be understood as follows. A uni- 
form phase gradient (induced here by a constant vector 
potential) generates a counterflow current in the super- 
fluid. One contribution to the current J2k v n r ( e fc + Vj ^-) 
originates from the change in the occupation of the 
quasi-particles due to the change of the energy disper- 
sion. To linear order in A the last expression becomes 
Jj- ^ k d ek n F ■ 2A from which the second term in expres- 
sion l[14p readily follows. A similar expression holds for 
the phase stiffness of a BCS superconductor, but with two 
differences. First, the velocity of a particle in a parabolic 
band is energy dependent and must be left under the 
integral. Second, the zero temperature phase stiffness 
naturally arises in a parabolic band due to the change in 
the velocity to v + A/m with m being the effective mass 
of the particle whereas in the graphene spontaneous- 
coherence case it is accumulated at the Dirac-model's 
cutoff wavevector, and conveniently captured[l[ by sub- 
tracting i D (A = 0). 

The influence of disorder on p s (T) follows mainly from 
the finite life time of definite-momentum quasi-particles. 
The impurities reduce the energy gap, thus reducing the 



energy gain due to pair condensation and correspond- 
ingly its sensitivity to phase gradients. The reduction of 
the KT temperature due to disorder is readily obtained 
from Eq. (fT2|) by setting the temperature to T KT and us- 
ing the relation p s (T KT )/T KT = 2/tt. The pairing po- 
tential A(T KT ) must be calculated self consistently using 
Eqs.jUl). 




FIG. 1: (Color online) KT temperature vs. disorder strength. 
Disorder is parameterized by a = 1/2t s T c o = S/T c q. The 
solid line in this figure corresponds to the clean system and the 
broken lines correspond to (top to bottom) a = 0.1,0.5,0.8. 
In graphene bilayers e F /T c o ~ lOfcfd in the strongly interact- 
ing regime. 

Figure[T]shows the reduction of the KT critical temper- 
ature of counterflow superfluids due to disorder. In this 
figure we plot the dependence of T KT on the Fermi energy 
for various strengths of disorder. All energies are scaled 
with T c q. The strength of disorder is parameterized by 
1/2t s T c o = S/T c q (~ nikpd/n in the strongly interact- 
ing limit). In the weakly interacting regime, e F /T c o ^> 1 
and the KT temperature is of order of T c . The effect of 
disorder on T KT can then be inferred from its effect on 
T c . In that regime A(Txt) is small hence the disorder 
will significantly reduce T KT . As the interaction increases 
so does T c o, however the ratio T KT /T c o decreases hence 
A(Tkt) is relatively large. Thus, in the strongly interact- 
ing regime the effect of disorder on the KT temperature 
is relatively weak. 

Discussion — Since the pioneering work of Abrikosov, 
Gorkov and Anderson it has been understood that both 
the critical temperature and the order parameter of s- 
wave superconductors are unaffected by a sufficiently low 
concentration of non-magnetic impurities (H}. This re- 
sult, known as Anderson's theorem [ljj], asserts that even 
in the presence of (non-magnetic) impurities two states 
related to one another by time reversal symmetry will 
pair and condense. On the other hand, as experimentally 
observed [3] and theoretically explained by Abrikosov- 
Gorkov theory, the values of of T c and A are suppressed 
by magnetic impurities that act to reduce the binding 
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energy of Cooper pairs. Later work demonstrated that 
Abrikosov-Gorkov theory applies to other circumstances 
in which some perturbation breaks time reversal symme- 
try, for example thin superconducting films subjected to 
a parallel magnetic field H, 17 1 or superconductors with 
an exchan ge f ield in the presence of strong spin-orbit in- 
teractions [18( . In this work we extended the range of ap- 
plicability of Abrikosov-Gorkov theory to dirty graphene 
bilayers. The extension to semi-conductor based bilayers 
is trivial. 

The analogy between the effects of disorder on super- 
conductivity in metals and on counterflow superfluidity 
in bilayer systems is evident if time reversal symmetry in 
the former is replaced by particle hole symmetry in the 
latter. Kramer's degeneracy in a non-magnetic supercon- 
ductor translates to particle-hole symmetry in the spec- 
trum of the bilayer system. Indeed, if we consider a ficti- 
tious scenario in which the disorder potentials of the two 
layers are perfectly anti-correlated, i.e. U\(r) = — [/^(r), 
then particle- hole symmetry is preserved, t s = — r D , the 
pair breaking parameter a vanishes, and both T c and A 
retain their clean system values. Generic disorder poten- 
tials break particle-hole symmetry and therefore provide 
a particle-hole pair breaking mechanism. The influence 
of disorder on the critical temperature is determined by 
the value of a — S/2T C . Disorder will have a negligible 
effect on T c as long as a <C 1 reinforcing the plausibility 
of condensation in the strongly interacting regime. 

Ideas similar to the ones used here may be applied to 
study disorder effects on exciton condensation in quan- 
tum Hall bilayers. In current experiments [n| l/2r s T c 
is of order of unity suggesting that incorporating disor- 
der is imperative for a correct estimate of the mean field 
critical temperature and of T KT . It is also interesting to 
consider the influence of magnetic impurities on a bilayer 
graphene system. Due to the spin and valley degeneracy 
of the graphene sheets the exciton pairing is SU(4) sym- 
metric. Put differently, four identical superfluids coexist 
in the system. A scattering event by a magnetic impu- 
rity will thus scatter an exciton from one superfluid to 
another. 
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